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Abstract 

The relativistic version of the J-matrix method for a scattering problem 
on the potential vanishing faster than the Coulomb one is formulated. As 
in the non-relativistic case it leads to a finite algebraic eigenvalue problem. 
The derived expression for the tangent of phase shift is simply related to the 
non-relativistic case formula and gives the latter as a limit case. It is due to 
the fact that the used basis set satisfies the "kinetic balance condition" . 
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I. INTRODUCTION 



The J-matrix method, introduced by Heller and Yamani [T|,^ and developed by Yamani and 
Fishman , is an example of an algebraic method in quantum scattering theory. Comparing 
with the algebraic variational theories the method has been shown to be free of the false 
resonances problem It has been used in construction of the Gauss quadrature of the 
continuum |^ (see also 0), the definition and analysis of a reproducing kernel in the context 
of Harris eigenvalues . Quite recently, it has been used in formulation of complex-scaling 
method and in development of the multi-channel Green's functions ^ by means of a 
complete basis. 

The crux of the method is representation of the Hamiltonian in a suitable non-orthogonal 
basis changing the differential scattering problem into the purely algebraic one. Thus far 
only the non-relativistic version of the method has been formulated. The aim of this paper is 
to develop the simple relativistic formulation of the method in its theoretical framework for 
potentials sufficiently regular at the origin and vanishing at infinity faster than the Coulomb 
one. 

II. NON-RELATIVISTIC J-MATRIX - RADIAL KINETIC ENERGY CASE 

First we briefly review the non-relativistic Jacobi matrix approach introduced in Refs. 
and extended in We recall only the case when the potential vanishing faster than 
the Coulomb one is involved, as we shall formulate the relativistic formalism for this kind of 
potential. The Coulomb case is much more complicated and it will be considered elswere. Let 
{0j^}^o be either Laguerre or Gaussian (Hermite) basis set. The explicit forms of both bases 
as well as some other formulas concerning non-relativistic problem (see Ref. |^) are collected 
in table I. Only the second basis, i.e., the Gaussian one, forms an orthogonal set, hence, 
in general, the notion of biorthonormality is needed. The set {0^}^o is biorthonormal to 
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{0^}^o with respect to the unitary scalar product if = /q°° 0^(Ar)0^(Ar)(ir = 5mn- 

Biorthonormal basis functions {(p^} ^^re also given in table I. The important feature of the 
sets {0^} is that the radial kinetic energy operator: 

^°"T = "2dH + ^^~T 

if expanded in any of them, takes the tridiagonal or Jacobi form: 

J^, = {<PU{Ho - ^V2)0L), Jmn + only for m = n, n ± 1. (2) 
In the above /c is a wave number related to the energy E and mass m of the projectile 

= (3) 

It must be stressed here that the matrix elements Jmn are functions of /c, i.e. J^n = Jmn{k). 
The regular solution S{k, r) of the equation 

{Ho-e/2)Sik,r)=0 (4) 

is simply proportional to the Riccati-Bessel function, satisfying S{k, r) ~ r'"*"^ as r ^ and 
S{k,r) ^-^^ sin(A;r — y). Using an expansion of S{k,r) in the basis {0^}, i.e. S{k,r) = 
Y^'^=o s''^(j)''j^{Xr) , one can write equation (|[) in the form 

oo 

E JmnSl = 0. (5) 
n=0 

As shown in Ref. using the explicit form of the matrix elements Jmn one can find the 
expansion coefficients in terms of Gegenbauer polynomials (see table I). Again we have 
^n — ^nW- J-matrix method to solve a scattering problem one introduces the 

second, cosine- like function C{k, r), which is required to satisfy C{k, r) ~ r'"*"^ as r — »• and 
C{k,r) '^^^ cos(fcr — y). It cannot be the second solution of the original, homogeneous 
problem as this solution, proportional to the Riccati-Neumann function, is singular at the 
origin. 



The required C{k,r) function has been found |T|-§[, in another way, namely by solving 
an inhomogeneous equation : 

{Ho-e/2)C{k,r)=0,{Xr), d = (6) 

with Sg being the first expansion coefficient of sine solution. Then the expansion coefficients 
of C(/c, r) satisfy the equation 

oo 

j^j^ = 0,. (7) 

n=0 

The corresponding coefficients = c^(fc) (see table I) have been also found by some 
differential technique . The calculated expansions S{k,r) = J2'^=o Sn't'ni^''^) C{k,r) = 
J2'^=o c^^nl-^^) have been used in an approximate solution of the original scattering problem 
on the radial potential V = V{r) vanishing faster then the Coulomb potential: 

{Ho + V- ^)iJE = 0. (8) 

Namely, the potential V has been replaced by a truncated potential operator 

= pI,VPn, (9) 

where P/v is the generalised projection operation: 

^^=El</'n)(0nl- (10) 
n=0 

The new potential operator can be written in the basis {</>^} as an x matrix with the 
matrix elements = (0jj|V^0^)- Then the exact solution ip'^ of the new problem: 

{H, + V^- = (11) 

has been expanded in the basis {0^} as 

N-l oo 

i^E^r) = E + E (^n + tan5^cl)0l (12) 

n=0 n=N 



to satisfy the boundary requirement ip^ij-) '^■^^ sin(A;r — y) + tan (5jv cos(A;r — y). The 
tan 5]^ is an approximation of the tangent of the sought phase shift 5 of the exact solution 
'ipE the problem (^. The left-hand side projection of (|Tl]) onto the basis {0^} gives then 
infinitely many equations depending on n. However all equations for n> N + 1 are satisfied 
automatically as coefficients sj^, cj^ satisfy the same recursion relation (|]) for any m > 0. 
The remaining finite set on equations involve + 1 unknowns tan^jv, {ani}n=Q- Those 
equations can be easily solved In particular, using the recursion relation for matrix 

elements Jnm the tangent can be calculated giving 

tan(5jv = — I ^ TTTTT (13) 

Ctv-I + gN-l,N-l\0}JN,N-lCi^ 

where gN-i,N-i{£-) = Y.n=d'^N-im/i^rn — £) with the matrix T diagonalising the finite- 
dimensional problem {T'^P'^{Ho + V — Y)P^)mn = i^n — £)Smn- Here the energy dependent 
quantity gN^i^N^i{S) can be viewed as the matrix element of the inverse of the truncated 
operator P^(//o + — if restricted to the A^-dimensional space where it does not 

vanish. The quantities £n [UHl (see also and references therein). 



III. RELATIVISTIC JACOBI-MATRIX PROBLEM 



Now we shall turn to the relativistic problem. Before the formulation of the method we 
shall find the relativistic counterparts of S{k, r) and C(fc, r) in some suitable basis. We shall 
also calculate the relativistic Jacobi matrix elements in this basis. For this purpose consider 
the free Dirac equation: 



/ 



[mc 



E)/ch -d/dr + n/r 



Fir 



^ d/dr + n/r {-mc^-E)/ch j \^G{r) j 

In the above the total energy E is related to the rest energy £ E = £ + mc^. Let 1{k) be 
the non-negative solution of the equation /(/ + 1) = n{n + 1), i.e. 1{k) = k, and 1{k) = —k — 1 



(14) 



for positive and negative k, respectively. We shall usually omit the symbol k in the notation 
throughout the text and write I only, remembering that the latter depends on k. Then 



equation (14) has two independent solutions. The first one, regular at the origin 



reg 



Gregir) 



J 



V 



(15) 



is constituted by the Riccati-Bessel functions with boundary behaviour 



.0 X- 



i+i 



and ji{x) 



(2/ + 1)!! 

The numbers e and k in (p!5f) are standard abbreviations 



sm x — 



(16) 



e = \i — k = 



E + mc^ ' 



ch 



The quantity k converges in the non-relativistic limit c oo to the number k 
second solution of (p!4D, irregular at zero is given by 



2m£ 



(17) 



The 



^irr y) 



±eni^i{kr) 



;i8) 



Here we have the Ricatti-Neumann functions with properties: 



,_o (2/-1)!! ^ _ , ™ , nl 
ni[x) — > ; and ni[x) — > — cos(x ) 



X 



(19) 



In both solutions (|T5|) and (|T8|) the upper and lower signs in the small components correspond 
to negative and positive k, respectively. From the above it can be immediately seen that 
the regular solution ^reg is the relativistic counterpart of non-relativistic function S{k,r). 
For the sake of consistency with the non-relativistic case, hereafter we shall denote the 
relativistic sine-like solution "^reg by "^s- To develop the Jacobi matrix analysis we have to 
introduce a suitable basis set. 
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A. The basis set 



in the Hilbert space L^(0, oo) ® on which the Dirac operator from (0) is defined. 
Let again be either the Laguerre or the Gaussian basis set and let ip^^^^Xr) = {n/r + 

d/dr)0^(Ar)[] Then the basis set defined for our purposes is 

/ . \ 



V 











(20) 



The above set depends on the positive reals number A which can be treated as a nonlinear 
variational parameter (see, for instance, [|I^. Note that the set (pOD satisfies the "kinetic 
balance condition". The latter condition is generally defined as a requirement that, if the 
funcitions {7^} are used to expand large component of solution of Dirac equation, then the 
basis {uji} used for expansion of small component should consist of linear combinations of 
functions {(^/r + d/dr)7i}. Use of such a basis is the simplest way |l^,14l to omit the 
problem of so called "finite basis set disease" (see HTBHTHl) in estimation of bound states 
of the atomic system. It seems that it would be also interesting in future to consider the 



relativistic J-matrix problem in the context of the relativistic Sturmian basis (see [19| and 
references therein) as it is known that the relativistic free particle Green function takes 
particularity simple form in this basis. 

The biorthonormal elements to the functions (pO]) obviously are <l>^^(r) = (0^(Ar), 0)"^, 
$~^(r) = (0, ^^(Ar))-^. As usual, we denote by /„ the element biorthonormal to /„. The 
elements 0^^"-' B are recalled in table I. Here we shall calculate the biorthonormal elements 



liK) 



One should keep in mind that here dependence on n is not only present via I coefficient, but via 
operator (^/r + d/dr). 
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It is easy to show integrating by parts, that biorthonormal elements {'0i(^)}^o should 
satisfy the equation 

/' 

Jo 

Hence it suffices only to solve the following inhomogeneous differential equation 

The resulting functions are given in Table II. They all belong to the space L^(0, oo). 
This fact is obvious apart from the case of negative k for the Gaussian set. This case needs 
more careful analysis as here it is not possible to give the functions by explicit formula. For 
all K < the functions {V'^} due to Gaussian set behave as r'+^ at the origin, and vanish 
at infinity not slower than r~^''^^^ as the limit of the occurring integral is finite. Then (V^^)^ 
behaves for r — > oo as r~'^(.''+^) and, as I is nonnegative, exists. Thus in both cases, 

when is either the Laguerre or the Gaussian basis set, all elements {V'^} biorthonormal 
to new functions ipl^ = {k/t + d/dr)4>l^ belong to L^(0, oo). Then obviously biorthonormal 
elements I»„k due to the relativistic case belong to the Hilbert space -L^(0, oo) (8) C^. Note 
that we do not need the explicit forms of biorthonormal functions in our considerations. 



dr 



^^(Ar) 



<Pl{Xr)d{r) 



(21) 



B. Expansions of relativistic sine and cosine solutions 

Now we are in the position to find the expansions of sine-like = ^^eg and cosine-like 
solutions. For the latter we demand to satisfy three requirements: 

(1) ^'c should have ^'j^r type asymptotic form, 

(2) should exhibit regular behaviour at the origin, 

(3) coefficients of '^c expansion should satisfy (apart from at most the few first ones) 
the same recurrence equations like the ones of "if reg = ^s- 
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Consider first the solution \l/(/(A;,r) = \^ Fi/{k,r), Gu{k,r) ] of tlie inliomogeneous 
equation of type (|I4D: 



[Ho - — )^C7 = <l>inh- 



(23) 



In tlie above the index U = S,C corresponds to sine-hne and cosine-hke solution. The 
inhomogeneity is chosen as ^inh = ^ok = ^u4>o, ^ and the coefficients Qu, U = S,C 
are fls = 0, = ~^/^o- 

Equation (^) can be also written as 

{n/r - d/dr)Gu - keFu = VLui^'^ 
We can introduce the relativistic counterpart of the Jacobi matrix : 



(24) 



j::i^{^U{'n,-E/cm:^:), s,s' = ±, m,n = 0,l,2, 



(25) 



The matrix elements of J can be expressed in an extremely simple form. To see this, let 
us define the 2x2 matrices Jmn defined by their matrix elements as {Jmn}ss' = Jmn- Then 
it can be easily seen that in the spinor basis the new matrix takes the particularly simple 
form: 



Cfmn 



(26) 



1^ (ClV'i) -WmWn) ) 

The explicit forms of the integrals constituting elements of the above matrix are given in 
table III. They are simply related to the non- relativistic J-matrix elements (|D(c.f. 0): 



(27) 



Now we shall predict the expansions of the two solutions in basis (^) in the following form 

/ „ \ 

^ U = S,C; u = s,c (2J 



E E = E ^nCk) 

S = ± 71=0 71=0 



V 



y 
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i.e. we predict that large components of sine-hke and cosine-hke solutions are given by the 
same expansion coefficients s^, cj^ as in the non-relativistic case, only taken in the modified 
point k and that the small components coefficients are only rescaled by e/k. 

It can be easily verified that ( pHf ) really solves the equation (|T^. Namely putting the 
above expansion into the equation and using the definition of matrix elements we get 
the infinite set of equations: 

oo 

E E ^r^nwl = ^f/fco5s,+, s = ±,m = 0,l,2,... (29) 

s'=± n=0 

as for any pair m, n fixed the second element in lower row of the matrix (|26|) is rescaled by — 
we obtain immediately that all equations (p9|) with a negative "index" s = "-" are satisfied 
trivially. Recalling the definition of / = /(/t) (cf. the remark following equation ([T^) ) after 



integration by parts one gets (t/^^K) = ((f + fMi^ + = (CK-jfr + 

Taking into account the form of the upper row of the matrix (p6D and the identity (|27|) we 
obtain immediately that equations (^9]) with the "index" s ="+" have the identical form 
with the sets of equations (|), if the latter are evaluated at k instead of k. Thus we 
have shown that the expansions (^8|) are in fact solutions of equation (p^. From the non- 
relativistic case we see that their large components have the desired behaviour at the origin 



and at infinity. Moreover, as the equations of the type (14) are coupled and the behaviour of 
the one component determines the behaviour of the other one. Hence both the components 
of the solutions \E'5, "^c have the asymptotic behaviour we need for purposes of our method, 
i.e., ^l/^ is simply the regular solution, behaves as ^jj-r at infinity and as "^reg at the 
origin. As the inhomogeneity involves only one biorthonormal element \E'(J"^ both functions 
satisfy the same set of equations apart from the first one (see formula (p6|)). 
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IV. POTENTIAL SCATTERING 



Now we shall consider the central problem of our paper which is the approximate solution 
within the jT-matrix formalism. Consider again the radial part of the scattering problem of 
a projectile on a target described by a sufficiently regular potential V = V{r) vanishing at 
infinity faster then the Coulomb potential. To solve the problem one has to find the solution 
of the following equation 

(^0 + ^ - ^) "^E = 0. (30) 

A solution of the above equation is required to satisfy the boundary condition 
E{k,r)^^^ s{k,r) + t^cik,r) where the tangent of the phase shift, t = tan5, is to 
be found. 

To develop the formalism of the relativistic J-matrix (we shall denote it by jT-matrix to 
distinguish from the non-relativistic case) we use the generalised projection operators: 

N-l 

Vn=Y.Y.\K.){K.I (31) 

s=± n=0 

and introduce the truncated potential: 

= '^N^'^N. (32) 

where Vj^ corresponds to hermitian conjugate of Vn- 

Now one can seek the exact solution of the equation with truncated potential: 

(Ho + V^-^)v^^(r) = 0. (33) 

Note that for any ip G L^(0, oo) (g) the function V^ip vanishes at infinity faster then 
^. Recall that we assumed that our original potentials vanish at infinity faster then ^. 
Thus although equation (|33D has not a standard Dirac equation form with the same scalar 
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potential in its large and small part, still its solution asymptotically satisfies free Dirac 
equation. Hence the solution \E'^ satisfies the boundary condition 



(34) 



ch 



where is an approximated tangent of phase shift. As the potential operator V 
we expect that for oo, converges to correct value t = tan 5. 

Now we shall find more details about the form of the solution The most general 
formula is 



(35) 



s=± m=0 

Consider the matrix representation of equation (^). Putting the expansion of the function 
in the basis {^1/^^} S^t the infinite set of equations: 



s'=± n=0 



N\ss' is' 
)mn nn 



0, 



±,m = 0,1,2, 



(36) 



It can be easily seen by the right-hand side projection of equations (^) onto the basis {^^k}- 

According to analysis following the formula ( pUD the expansion coefficients {d^^} of 
must satisfy: (a) for the large component Z^^o Jmn{}i)dnK = 0, m > N with elements Jmn(-) 
given by the non-relativistic formula, (b) for the small component d^^ = Moreover 
we impose the additional condition (c) ~ S{k,r) + iiyC{k,r) (see condition (|34D ). This 
gives us, together with the condition (b), the following required form of the sought solution 
of equation (||) (c.f. [|): 



N-l 

E 

m=0 



+ E 



m=N 



Km /rn 



V 



+ tNC 



(37) 



where the abbreviations has been used according to (p8|). After adding and sub- 

T 



tracting the term Ylm=o 



Km/T^n^ \" Km 



to the left hand side of the 
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above equation it is straightforward to see that the above function satisfies the asymptotic 



condition (p^). 

Let us turn back to equations In general, in analogy to the non-relativistic case, 

they can be schematically represented as follows: 



/ 



XXX 
XXX 
XXX 
XXX 







XXX 
XXX 

X X X X X X 

X X X X X 

X X X X X XXX 

xxxxxxxx 

X X X X X X X 
X X X X X X X 

xxxxxxxx 

XXX X X X X X 

XXX X X X X X 

XXX X X X X X X 
XXX XXX 

XXX XXX 



\ ( 



V ■ ■ 



• • 7 



iV+l,i 
s+ 

i+ 



4,1 



N-l.l 
*JV+l,i 



V 







V7 



From the construction of the required form (pTf ) we see that all the equations for m > N 
are satisfied automatically. Thus one has to solve the remaining 2N + 2 equations with the 
unknowns i^, rfg^, df^^, d^^, di,^, Note that here the number of equations 
is greater then the number of sought quantities {2N + 1), so in general the set of equations 
of such a form can have no solution. But in our particular case the solution certainly exists 
as the general theory of differential equations assures the existence of and, according to 
the previous analysis, ( p7| ) represents the most general required form of 

Using equations (|36D one obtains the following form of the remaining equations: 



13 



/ 7++ 







JV,JV-l^JV-l •^JV,JV-1 



JV,]V-1 



JV-l,JV-2 



J, 



JV-1,]V-1 



JV-2,JV-2 



JV-2,]V-1 



(^+V^)oV2 
(^+V^)0,^-2 



- 7+" c" \ 

-'jV,JV-l'-JV-l 
•-'jV-l.iV^JV 








7"+ c+ 



JV,JV-1^JV-1 



JV-2,JV-1 



JV-1,JV-1 

7"+ 



^ N- 


4- 

-2,JV- 




f 

-1,JV- 




( ) 




4-1 




"jV-2 




4 








'^JV-2 








V / 





JV-2, JV-2 (•^ + "'^"'^)jV-2,JV-l 
JV-l.JV-2 -"JV-I.JV-I 



J: 



JV,JV-1 



•~^n-i,n'^n 

~^ /V /V — 



-l/ 



- 7++ s+ - 7+" 







JV-l.JV^JV 



JV-I.JVJV 



JV,JV-1 JV-1 



JV,JV-1 JV- 



Keeping in mind that the inner 2N x 2N matrix {J^+V 



N\ss' 



±,m,n = 0, 1, 



is Hermitian and real, hence symmetric, and recalhng the definition of we can solve the 
above equations by some orthogonal matrix F (cf. the non-relativistic case): 



2N X 2N matrix G{E) with elements defined as 

N-l ysp ys'p 



(3J 



p=± i=o - ^ 



(39) 



is an inverse of the 2N x 2N matrix representation of the truncated operator V\q{7iQ + ^ ~ 
^)Vn- It can be viewed as the approximation of the relativistic Green function in the basis 
(|20|) . The numbers are the relativistic counterparts of the Harris eigenvalues They 
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represent a finite approximation of the spectrum of the relativistic Hamiltonian Tip + ^• 
In particular, they include positive approximations of the first energy levels due to the 
potential V and the negative pseudo-energies due to the continuous spectrum. As our 
basis ([20| ) satisfies the kinetic balance condition [0,0] there is a hope that Ef satisfy the 
generalised form of the Hylleraas-Undheim theorem (see, for instance [p!5| , |20| and references 
therein). It means, in particular, that (i) N positive values among set {Ef} approximate 
the exact eigenenergies Ho + ^ from the above and that (ii) the remaining N eigenvalues 
have values below —mc^. 

We can introduce now the {2N + 2) x (2N + 2) block-diagonal matrix : 

r(2Af+2)x(2Ar+2) = diag{l, T2NX2N, 1) (40) 

and act with it on the left-hand side of the above set of 2A^ -|- 2 equations. Using the fact 
that the matrix i7iv,Af-i given by (p6|) is nonsingular the set of (2A^ -|- 2) x (2A^ -|- 2) the 
equations can be solved with respect to the approximate tangent of phase shift. Using the 
properties of the coefficients of the matrix J' one can derive the tangent of the appoximated 
phase shift in the form similar to the non-relativistic formula: 

^ s^N-iCk) + (2e/fc)g+li^^_i(E) J;v,^-i(fc)4(fc) 

4_i(^) + (2e/^)6?+li,^_i(i5;)J;v,iV-i(^)4(^)' 

Note that in the above the Jn Ar_i stands for the non-relativistic J-matrix element (see 



(|27|) ). The fact that we have 2A^ -|- 2 equations and 2A^ + 1 unknows results in second, very 
similar formula for t^r with (A;/e)^^l^ ^_]^(£') instead of G^'ti n^i{E). From the previous 
analysis we know that both equations must give the same which means that one has 
GnU^n-AE) = {e/~k)g^t.^.,{E). 
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V. DISCUSSION 



Comparing equation (^iD with the non-relativistic formula (0) one can see that apart 
from the quantity {2e/k)Q^'ti all elements of the expression for tangent of the phase 

shift have the same form as in (|13D, they are only evaluated in relativistic wave number k. 

Now let us note that for any N the above formula for tangent shift converges to the 
non-relativistic limit as the speed of light c approaches infinity. Indeed, the used basis ( pOD 



ensures see 



1^,|T5[) that in the limit of infinite c the large component satisfies the correct 
Schrodinger equation ( pT]) with the wave number k = limc^oo k. This means that the related 
tangent of the phase shift must also satisfy a correct limit, i.e. 

lim In = tjy. (42) 

c— >oo 

From the above we get immediately \imc^oo{'^^/k)Q^'tii^_i{E) = gN-i,N-i{^)- Moreover 
{2e/ k)Q^'ti N^i{E) plays the analogous role as gN-i,N~i{^)- In fact, the matrices Q{E) and 
g{S) can be viewed as the finite approximations of the Green functions of the relativistic and 
non-relativistic Hamiltonians with the potential V, respectively. The form of the factor 2e/k 
is simply connected with the normalisations of the Green functions in both cases. It can be 
seen from the simple analysis of the set of second order equations derived in a standard way 
from the Dirac equation. 

From the practical point of view, the convergence can be improved with the help of 
additional parameter A. As we mentioned before, the latter can be treated as an additional 
variational parameter. In particular its optimal value will depend on the range of the 
potential. It can be simply seen that potentials of long range should be treated with small 
A while potentials with support located close to the origin will require large values of the 
parameter. 

In conclusion, we have provided the relativistic version of Jacobi matrix method for well 
defined class of potentials. The usage of the basis satisfying the "kinetic balance condition" 
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allowed for a simple formulation of the method. In particular, the derived expression for the 
tangent of the phase shift is similar to its non-relativistic counterpart and reproduces the 
latter as a correct non-relativistic limit. 
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on kinetic balance condition and P. Syty for remarks on the manuscript. The work is 
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TABLES 

TABLE L Elements of expansions of sine- and cosine-like solutions in the Laguerre and the 
Gaussian basis sets The L^n^ and ci"^ are the Laguerre and the Gegenbauer polynomials, 
respectively while 2-^1 and \F\ are the Gauss and the Kummer (confluent) hyper geometric functions 
(see [^) respectively; A > is a scaling parameter. 



Quantity Laguerre set Gaussian set 

(/-J, (Ar)'+i exp(-Ar/2) L^n^^^ (Ar) (Ar)'+i exp(- AVV2) L^n^^^ (AV^) 



An.+2«+l (^'') '^n A^r(n+2i+3/2) '^n 



4 (^^^ exp(-,V2)Lr^/^^ (,^) 

0Fr(n+2«+2)(sine)' r(n+i+f) ^^'^ ' I > > 

X2Fi(-n - 2/ - 1, n + 1; i - /; sin2(e/2)), XiFi(-n - / - i, i - /, t?^), 



TABLE IL Biorthonormal elements '0n due to the small component. 

Quantity Laguerre set Gaussian set 

for K > __^(Ar)^exp(-Ar/2) exp(-A2rV2)(Ar)' 

¥n for K < --^exp(-Ar/2)ELo(-2)^+^ - aT (n+lz+3/2) (^")'^'^'^ 

X Eio (g^^(Ar)^'-^-n-l)^4^r^^^(Ar) x J^)^ t'+V^ exp(-tV2)Lr^/^) (t)dt 
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TABLE III. The overlap integrals proportional to the elements of j7-matrix. 



Integral Laguerre set Gaussian set 

-{n + 2n + 2,)5m,n+i] 

{^LWn) ^^^^^^[2(n + 21 + 2)(2A - A^r(n+j+3/2) ^^^^ ^ ^ ^ 3/2)<5^„+ 

n(5„i,n-l + (n + 2?7. + 3)5r„,n+l] f^i^m.n-l + {n + I + 3/2)5^,„+i] 
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